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NOL HYPERSONIC TUNNEL NO. 4 RESULTS VII: 
EXPERIMENTAL INVESTIGATION OF TURBULENT 
BOUNDARY LAYERS IN HYPERSONIC FLOW 


Prepared by: 

R. Kenneth Lobb, Eva M. Winkler, and Jerome Persh 


ABSTRACT: Naturally turbulent boundary layers on the wall 
of a wedge-type water-cooled nozzle in the NOL 12 x 12 cm 
Hypersonic Tunnel No. 4 have been investigated at Mach 
numbers of 5.0 to 8.2 with and without steady state heat 
transfer to the surface. The Reynolds number based on 
boundary layer momentum thickness v/as varied from 5,000 to 
13,000. Measurements of Pitot and static pressures, total 
and wall temperatures, and rates of heat transfer made it 
possible to compute velocity profiles, temperature profiles, 
and boundary layer parameters without resorting to any 
assumptions. The turbulent portion of the boundary layer 
velocity profile was found to differ from the incompressille 
flow logarithmic law by an amount that depends on the heat 
transfer and Mach number. The data for the outer turbulent 
portion for any one Mach number fall on a single curve if 
plotted in a particular nondimensional coordinate system 
which is based on local properties in the boundary layer. 

The velocity profile in the laminar sublayer is linear. The 
thermal sublayer was found in all cases to be larger than 
the velocity sublayer. Local skin friction coefficients 
for zero heat transfer as calculated from the velocity gradi' 
ents at the wall are consistent with the results of other 
experimenters at lower Mach numbers for the same Reynolds 
number. Furthermore it is demonstrated that the Reynolds 
analogy between skin friction and heat transfer is valid at 
Mach numbers up to 8.2. 
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This is the seventh NAVORD Report on an investigation 
carried out in the continuous NOL 12 x 12 cm Hypersonic 
Tunnel No. 4. The titles of the previous NAVORD’s describing 
results from the tunnel are: 

I. Air Liquefaction. NAVORD 1742, 4 January 1951 

II. Diffuser Investigation. NAVORD 2376, 5 Hay 1952 

III. Diffuser Investigation with Models and Support. 
NAVORD 2435, 1 July 1952 

IV. High Supply Temperature Measurement and Control. 
NAVORD 2574, 8 October 1952 

V. Experimental and Theoretical Investigation of a 
Cooled Hypersonic Wedge Nozzle. NAVORD 2701, 

13 April 1953 

VI. Experimental and Theoretical Investigation of the 
Boundary Layer and Heat Transfer Characteristics 
of a Cooled Hypersonic Wedge Nozzle at a Mach 
Number of 5.5. NAVORD 3757, 8 July 1954 

The present NAVORD (Results VII) presents experimental 
investigations and discussions of turbulent boundary layers 
in hypersonic flow. 

Knowledge of the related effect of skin friction and atro- 
dynamic heating at high Mach numbers is needed by the 
designers of hypersonic vehicles because, for slenc’er 
missiles, skin friction represents the major part of the 
drag and extreme surface temperatures produced by frictional 
heating may very well be a limiting design factor. 

A portion of the results contained in this NAVORD were 
presented at the 22nd Annual Meeting of the Institute of 
Aeronautical Sciences and also at the Bureau of Ordnance 
Committee on Aeroballistics Symposium in October 1954. The 
present report contains additional results and more detailed 
analysis of the data as well as a complete tabulation of the 
experimental results. 

This work was jointly sponsored by the U. S. Naval Bureau of 
Ordnance an«’ the U. S. Air Force. It was carried out under 
Tasks N0L-M9a-108-l-54, NOL-M9a-133-l-55. and 
N0L-M9a-133-5-55. 

The authors are indebted to Dr. R. E. Wilson for many stimu¬ 
lating discussions during the course of the investigations. 
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The numerical evaluation of the boundary layer surveys was 
done on a card-programmed calculator of the Applied Mathe¬ 
matics Division. The preparation of the data for the card- 
programmed calculator was greatly assisted by Dr. E. K. Blum. 
The cooperation of Messrs. L. L. Liceini and R. Carreu, Jr., 
who participated during the tests, is acknowledged. 


JOHN T. HAYWARD 
Captain, USN 
Commander 


II. H. KURZWEG, Chief 
Aeroballistic Research Department 
By direction 
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SYMBOLS 

Icon 1 skin 
conditions 

incompressible local skin friction coefficient for 
zero heat transfer based on froe-streara conditions 

specific heat at constant pressure 

boundary layer shape parameter ^*/6 

thermal conductivity 

Mac.. umber 

Nusselt . -or for stagnation temperature probe 
Prandtl numbei 
stagnation pressure 

Pitot pressure 

static pressure 
Gas constant 

Reynolds number based on free-stream conditions 
recovery factor 
Stanton number (equation (3)) 
local stagnation temperature 

stagnation temperature as measured by a stagnation 
temperature probe 

local static temperature 

velocity 

velocity component perpendicular to wall 
friction velocity (T^/p) 2 


friction coefficient based on free-streas 


go u o 
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u + velocity parameter u/uf (p based on wall properties) 

u* velocity parameter u/u^- (p based on local properties) 

y distance perpendicular to wall 

y* wall distance parameter y /v* ( p and P based on 
wall properties) 

y* wall distance parameter y Uy /** (p and based on 
local properties) 

f ratio of specific heats 

<f total boundary layer thickness 

q * boundary layer displacement thickness 

<r 

9 momentum thickness / — P — —fl - ^—]dy 

-f Jp00 U 00 L 



/* 

9 

r 


viscosity 

kinematic viscosity 

density 

shear stress 


Subscripts 

e equilibrium wall temperature for zero heat transfer 
g physical properties of gas 

L edge of laminar sublayer 

p temperature probe 

L physical properties of thermocouple wire 

T temperature profile 

u velocity profile 

w values based on wall conditions 

9 values based on momentum thickness 
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values bar'd on free-stream conditions outside 
boundary layer 

conditions behind shock in front of probe 
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.'OL HYPERSONIC TUNNEL NO. 4 RESULTS VII: 
EXPERIMENTAL INVESTIGATION OF TURBULENT 
BOUNDARY LAYERS IN HYPERSONIC FLOW 


INTRODUCTION 

1. Boundary layer investigations at hypersonic speeds are 
of immediate practical interest since friction drag and 
heat transfer data are needed by the designers of hypersonic 
missiles. 

2. From the theoretical viewpoint, various treatments of the 
problem are available for the case of laminar boundary layer 
flows as well as for turbulent boundary layers. The experi¬ 
mental data are, however, limited, and in general cover the 
Mach number range up to about 5. The lack of experimental 
confirmation for laminar boundary layer theory is not 
regarded as serious since the theory is considered to be more 
or less exact. Most of the analyses of compressible turbu¬ 
lent boundary layer flows, however, arc based on experimental 
results obtained in wholly incompressible flow. Tne uncer¬ 
tainties in these analyses have led to large discrepancies 

in the prediction of skin friction and heat transfer. Even 
those theories which agree with each other and with the exist¬ 
ing experimental data within 5 percent up to Mach numbers of 
5 differ greatly at hypersonic Mach numbers {reference (a)). 

3. The purpose of the present investigation was to exten.o 
the Mach number range of available data and also to attempt 
to provide a deeper insight into the characteristics of a 
turbulent boundary layer m compressible flows. 


Experimental Equipment and Techniques 

4. The boundary layer surveys have been conducted in the 
NOl 12 x 12 cm Hypersonic Tunnel No. 4, which is described in 
references (b), (c), and (d). This tunnel operates contin¬ 
uously m the Mach number range from 5 to 10 free of air 
condensation effects. Supply temperatures from 300 to 800° K 
and supply pressures from 1 to »5 atmospheres are available. 

A two-dimensional water-cooled wedge nozzle expands the air 
and can be adjusted for any Mach number simply by changing 
the throat opening area. 

5. The investigations were made for a ran ,c of Mach numbers 
varying from 5 to about d. For each Mach number several 
surveys were mauc diffenn.. *n the Reynolds number and the 
rate of heat transfer to the wall. The surveys were made 

at the center of one nozzle wall, approximately four inches 
upstream of th<- nozzle exit plane. Figure 1. Measurements 
at this pos.l.oa are unaffected by the junction of the nozzle 
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c,io' and the first diffuser plate (reference (e)). The rela¬ 
tively thick (approximately 25 mm) turbulent boundary layer 
on the nozzle wall minimizes probe positioning errors and 
thus facilitates accurate determination of the profile shapes. 
Such a boundary layer, however, is subjected to a slight 
frec-strcam pressure gradient in the neighborhood of the 
survey station because the flow is radial. (The corresponding 
Mach number rise is about 3 percent per tunnel caliber at the 
survey plane for a Mach number of 5 and decreases with 
increasing Mach number.) it is felt that the effect of the 
pressure gradient on the profiles is small. The wall tempera¬ 
ture, i the other hand, is maintained practically constant 
near room tempera Lure over the entire length of the nozzle, 
except at the throat where the surface temperature approaches 
the recovery temperature. In all tests the transition from 
laminar to turbulent boundary layer occurred slightly down¬ 
stream of the nozzle throat*. The exact effect of the history 
on the local boundary layer characteristics is, of course, 
not fully known and the results must be interpreted with this 
fact in mind. 


6. Reynolds numbers and the rate of heat transfer to the 
wall are controlled by the supply pressure and the supply 
temperature. In all tests the Reynolds numbers based on 
boundary layer thickness are roughly two orders of magni¬ 
tude greater than needed to make slip flow effects at the 
surface negligible (reference (f)). Since the wall tempera¬ 
ture is always maintained near room temperature, the lowest 
rate of heat transfer to the water-cooled nozzle that could 
;.<■ realized corresponds to a supply temperature which is just 
high enough to avoid air condensation in the test section. 

7. For each survey. Pitot and static pressures, stagnation 
and wall temperatures, and the wail temperature gradient 
perpendicular to the nozzle surface are recorded. Static 
pressures are measured by a 0.G4 mm diameter orifice in the 
wall and by a static probe in the free stream just outside 
the edge of the boundary layer. The static probe is an 8° 
cone cylinder with orifices located 17 diameters aft of the 
shoulder, Figure 2. It is mounted from the side wall with 
its axis parallel to the flow. Since agreement between 
wall- and free-stream static pressures is within 1 percent, 

a constant static pressure has been assumed to exist through 
the boundary layer. The measurements arc made with silicone 
oil manometers of±2 microns measuring accuracy. 


♦Surface probe tests made subsequent to the results reported 
jr. reference (e) indicated that the boundary layer profile 
reported In this reference was close to the transition region. 
It was definitely established, from the results of the sur¬ 
face probe tests, that all boundary layer profiles presented 
herein were measured in the fully turbulent .region. 
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8. The Pitot pressure is surveyed from wall to free stream 
with a flattened hypodermic tube of 0.125 mm half height, 
Figure 2. The opening is large enough to avoid errors due 
to slip flow effects on the pressure measurements in the 
region of low Reynolds number in the boundary layer near the 
wall (reference (g)). The probe is mounted in a micro¬ 
traverse mechanism which also accommodates the connection 

to the pressure gauge. The position of the centerline of 
the probe is measured relative to the position of electrical 
contact between probe and wall with an accuracy of ±0.025 mm. 
Impact pressures above 20 ram Hg are measured with a precision 
mercury manometer of±0.1 mm measuring accuracy. For the 
lower pressures a silicone oil manometer is used. 

9. A stagnation temperature probe with a single platinum- 
coated silica shield is used for the temperature surveys 
of the boundary layer, Figure 2 (reference (h)). For the 
rreasurements close to the wall, a flattened probe with a 
half height of 0.48 mm is used. The temperature recovery 
factor of tho probes reaches a value of 0.998 for large 
Reynolds numbers. The variation of the temperature recovery 
factor with flow parameters for each probe is described by 

a single calibration curve by relating the calibration data 
to the flow conditions inside the probe. Reference to this 
curve makes it possible to determine accurately total and 
static temperatures throughout the boundary layer from 
measured temperatures and pressures. The e.m.f. output of 
the temperature probe is recorded on Brown temperature 
recorders. The temperature measurements are accurate to 
± 0 . 2 ° 0 . 


10. Local values of the heat transfer to the nozzle wall and 
nozzle surface temperaturos at the boundary layer survey 
station are obtained from temperature measurements in tho 
nozzle wall. These measurements are made with four thermo¬ 
couples imbedded in the nozzle wall at various distances 
from the surface (reference (d)). Depending upon the 
operational conditions of the tunnel and the rate of coolant 
flow to the nozzle, equilibrium readings are reached after 
10 to 20 minutes of tunnel operation. They show that the 
temperature drops linearly from the nozzle surface. Previous 
investigations indicate that longitudinal and lateral tempera¬ 
ture gradients at the nozzle surface can be neglected 
(reference (d)). The wall temperatures are measured to within 
±0.01° C with a Leeds and Northrup K-2 type potentiometer. 


Data Reduction 

11. The experimental data needed to evaluate the boundary 
layer profiles are: 

a. Pitot pressure, p^, vs. distance from wall 
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b. Static pressure, p. taken constant 

c. Measured total temperature, T if vs. distance from 
wall 

d. Variation of probe temperature recovery factor with 
flow parameters 

Since p^ and Tj are not always measured at identical 

positions from the wall, both sets of data are taken at 
sufficiently close intervals so that the temperatures corre¬ 
sponding to impact pressure readings at a particular position 
can be interpolated from a curve faired through the tempera¬ 
ture data. 

12. The numerical evaluation of the data is done on the HOL 
Card-Programmed Calculator. The Rayleigh formula is used to 
compute the Mach number profile from p£ and p (reference 

(i)). To evaluate the stagnation temperature, Tq, from the 
measured T A values requires reference to tho calibration 

curve of the probe (reference (h)). This curve gives the 
recovery factor of the probe in terms of the measured total 
temperature and the pressure, p 2> behind the bow shock in 

front of the probe 


r 


P 


7/4 

f(P 2 /RTj ) 


( 1 ) 


13. To calculate the friction velocity, u^, for the u + , 
y+ representation (reference (j)) of the boundary layer 
velocity profiles, the wall shear stress, tT w , is calculated 

either from the slope of the velocity profile in the laminar 
sublayer 



( 2 ) 


or from tho temperature measurements in the nozzle wall, 
using the following equations: 


St 


uo 


( AT/Ay)k _ 

<T 0 " C p f i CT> U oo 


(3) 
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and the Reynolds analogy (reference (k)) 


-r 


St 


CD 


2/3 2 

** />CO 


(4) 


Values of 4T/Ay are obtained from the measured wall 
temperature gradient and values of T« were calculated using 
1/3 

r = Pr with Pr - 0.72. 


RESULTS 


Roundary Layer Profiles 

14. Naturally turbulent boundary layers were surveyed at 
free-stream Mach numbers of 5.0 to 8.2. These surveys were 
made at different Reynolds numbers Reg and various rates 

of heat transfer to the nozzle wall, as listed in the 
following table: 


00 

Re e 

(T e - T w )/T e 

4.93 

5,350 

0 

5.01 

6,480 

0.223 

5.03 

7,950 

0.373 

5.06 

7,370 

0.420 

5.75 

11,600 

0.108 

5.79 

12,400 

0.239 

5.82 

11,400 

0.379 

6.83 

8,550 

0.325 

6.83 

12,640 

0.443 

6.78 

8,400 

0.437 

6.78 

7,900 

0.500 

7.67 

8,130 

0.487 

8.18 

9,540 

0.496 


A tabulation of the complete data is given in Tabic I 
(page 22) and Appendix A. Typical plots of the basic data 
are shown in Figures 3 through 6. The data points close to 
the wall that wore influenced by the presence of the wall 
were omitted from all figures and tables. 
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15. For most of the curves, a dimensional distance has been 
selected for the aoscissa in order to show more clearly the 
physical differences between the profiles at different heat 
transfer rates. For the same reason some of the curves and 
data points have ueen omitted from Figures 3 through 5. Shape 
differences are pronounced only in the temperature profiles 
of the houndary layer, Figures 4 and 5. The static tempera¬ 
ture profiles, Figure 5, obtained with heat transfer have a 
temperature maximum close to the wall. The temperature at 
the maximum point or those still closer to the wall could 
not be evaluated from measured air temperatures because of 
the physical size of the temperature probes. In the cases 
where heat transfer data have been measured in the nozzle 
wall (see Table I), the slope of the temperature curve 
immediately at the wall has been deduced from the tempera¬ 
ture gradient in the nozzle wall 


VV«&JL X 


(5) 


(Values of the thermal conductivity of air have been taken 
from reference (a)). The curve through the measured data 
and the computed slope have been joined to give a smooth 
shape to the static temperature curves. (Since only the 
square root of the temperature enters into the computation 
of the velocity, the velocity profile close to the wall was 
found to be insensitive to errors in T due to incorrect 
interpolation between the last measured air temperature and 
the wall temperature.) 

16. The velocity profiles at any one Mach number, Figure 6, 
are similar in the turbulent outer part. A sharp change in 
slope, near the wall, roughly specifies the extent of the 
laminar sublayer. Within this sublayer the velocity varies 
linearly with distance; the slope, however, increases with 
increasing heat transfer. 

37. It should be noted here that velocities deduced from 
Pitot measurements in the laminar sublayer may be slightly 
high due to the effect of velocity fluctuations on the Fitot 
pressure readings cl ose to the wall. Since the fluctuating 
velocity component, uis always positive, the Pitot tube 
measures total pressures that are always somewhat greater 

than the effective total pressure by the amount lp u'^. 

~~2 2 it 

In Figure 7 the variation of (u' /u ) 2 with nondimensional 
distance, y*, from wall is shown for the incompressible data 
of references (m), (n), and (o). Although there is consid¬ 
erable scatter, the curve was drawn by giving the greatest 
consideration to the trend indicated by the larger bulk of 
data. On the assumption that the incompressible results of 
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Figure 7 are applicable to the present data, the laminar 
sublayer velocity profile data can be adjusted by the method 
given in reference (m). This adjustment indicates that the 
measured velocities at the edge of the laminar sublayer may 
be higher than the true velocity by a maximum amount of 
3 percent 


Nondimensional Velocity Profile Representation 

18. Turbulent boundary layer velocity profiles are very 
often presented in nondimensional coordinates which may be 
computed either on the basis of wall properties, u * and 
y"**, or local properties at each point in the boundary layer, 
u* and y* (see list of symbols). Figures 8 to 11 show 

all the measured velocity profiles in the u - **, y* coordinate 
system. A pronounced upward shift of the turbulent portion 
of the curves with increasing heat transfer is apparent for 
each Mach number. A similar shift has been observed by 
Deissler (reference (p)) in his experiments with subsonic 
compressible turbulent boundary layers in pipes. The position 
of the three upper curve* of Figure 10 relative to each other 
indicates that a Reynolds number effect is probably super¬ 
imposed. According to the respective values of the heat 
transfer parameters for these surveys, a different spacing 
in the turbulent portion would be expected. 

19. In addition to this Reynolds number effect the data in 
the outer turbulent portion of the profiles (figures 8 to il) 
apparently disperse with Mach number and heat transfer 
parameter, or some combination thereof. In an effort to 
determine a parameter which would eliminate the dispersing 
effect of at least one of the quantities mentioned above, 
the values of y» and y* were evaluated for each point 

in the boundary layer for each of the measured profiles. A 
plot of these results indicated a tendency toward a single 
curve for each Mach number, with the dispersing effect of 
the different heat transfer rates removed, but not that of 
the Reynolds number. It was found that this dispersion of 
the data could be removed by dividing the y* coordinate 
by the respective values of the shape parameter (H). The 
data for any one Mach number then converge to a single 
curve, (Figure 12). With increasing Mach number the curves 
fan out in an upward direction. That this displacement is 
systematic is demonstrated in Figure 13 which shows a cross 
plot of the data of Figure 12. The naturally turbulent 
boundary layer data of reference (t), as well as the M « 0 
results are included in Figure 12. 

20. The data of Figures 0 to 11 may be used to examine some 
of the fundamental assumptions of several of the theoretical 
treatments of compressible turbulent boundary layers. It is 
implied in these analyses (references (q), (r), (s) that the 
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edge of the laminar sublayer occurs at a fixed value of 
u* - y*, and that the slope of the turbulent portion of the 
u+, y* curve is constant, regardless of Mach number or heat 
transfer conditions. It should be noted that the results 
of Figures 8 to 11 show that these assumptions are not 
generally valid. 


Velocity and Thermal Sublayers 

21. It is apparent from the curves of Figures 8 through 11 
that the u*, y + velocity profile representation only roughly 
specifies the extent of the laminar sublayer. The edge of 
the velocity sublayer can be accurately determined from a 
plot of the coordinates u/uoo against y /<f on logarithmic 
paper. On such a plot it is found that the outer turbulent 
position of the boundary layer may be well fitted with a 
straight line, the slope of which is the exponent in the 
power profile representation of the turbulent boundary layer 
profiles 



<y/«n 1/r * 


( 6 ) 


In the laminar sublayer, a straight line of unit slope may 
be drawn quite accurately through the data. The point of 
intersection of these two straight lines is defined as the 
edge of the laminar sublayer. Figure 1*1 demonstrates the 
results of this procedure applied to the data of several 
typical profiles. The intersection points indicated are 
fairly insensitive to the slope of the line faired through 
the experimental data in the turbulent portion of the veloc¬ 
ity profile,and therefore the scatter of the experimental 
data does not seriously affect the determination of this 
point. A similar procedure may be used to determine the 
edge of the thermal sublayer. The thicknesses of the veloc¬ 
ity and thermal sublayers have been determined for all of 
the measured profiles and the values of u,/u_ , (ef./ef) , 

(c f /<f) , and n are tabulated in Tabl® I (page 22). 

L T 

22. A comparison between the values of (cT,/<f) and 

Li U 

( cf L /cO T , Table I, shows that the thermal sublayer is m 

all cases thicker than the velocity sublayer. According to 
the analysis of reference (v) the relative thicknesses of 
the thermal and velocity sublayers are related to the curbu- 
lent Prandtl num/cr. For the present case of thermal sub¬ 
layers thicker than the velocity sublayer Reichardt’s tneory 
predicts a turbulent Prandtl number larger than the molecular 
Piandtl number. Preliminary evaluation of this quantity fov 
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the zero heat transfer case of M - 4.93 and a molecular 
Prnndtl number of 0.72 using the equations given by Rubesiu 
(reference (w>) yielded a value for the turbulent Irandtl 
number of the order of 0.9. 


Discussion of Skin Friction Data 

23. Because many of the theoretical treatments of compress¬ 
ible turbulent boundary layers involve properties of the 
laminar sublayer, emphasis was placed on determining the 
velocity and temperature distributions in this region. 
Furthermore these data facilitate the accurate determination 
of surface shear stress from the velocity profile slope near 
the wall. 

24. For most of the measured profiles, the surface shear 
stress was calculated from the slope of the velocity profile 
in the laminar sublayer and the heat transfer data using 
Reynolds analogy, as described previously. Values of c f 
and St^ obtained for all the profiles arc tabulated in 
Table I. In general the values of c f determined from the 

velocity profile data are up to 5 percent larger than those 
determined from the heat transfer measurements. A possible 
explanation for this discrepancy way be given by considering 
the previously mentioned effect of the velocity fluctuations 
on the Pitot pressure data. Assuming that the incompressible 
results of Figure 7 can bo applied, the laminar sublayer 
velocity profile data were corrected by the method indicated 
in reference (m>. Values of wall shear stress were then 
computed from the adjusted velocity profiles. This procedure 
brings the shear stress values from the heat transfer data 
and velocity profile data in closer agreement (see Table I). 
The two values are well within the experimental accuracy of 
the data. 

25. The good agreement between skin friction coefficients 
obtained from the velocity profile data and those determined 
from the measured heat transfer, using Reynolds analogy 
demonstrates the applicability of Reynolds analogy at hyper¬ 
sonic Mach numbers. That the Reynolds analogy is applicable 
for supersonic flow up to a Mach number of about 3.2 has been 
shown by Seiff (reference (x)). The present data, plotted 
together with those of reference (x), show that at hypersonic 
Mach numbers the Reynolds analogy maintains the same form 
found to be accurate for lower supersonic speeds. The present 
data are especially convincing evidence because the Stanton 
number and the skin friction coefficient were ''etormmed 
simultaneously by two independent experimental methods whereas 
the results of reference (.\) were calculated by "interpolating 
the skin fnci ion data l«> equal couth i ioim of Mach number, 
wall to free stream temperature ratio, and Reynolds number.■ 
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26. In the figures which follow, the experimental values 
of c f used are those determined by the velocity-profile 

slope technique. These values of c^ were used, not in 

preference to those obtained from the heat transfer measure¬ 
ments , but only because there arc some cases for which heat 
transfer measurements are not available and it was desired 
to preserve consistency in the plotted results. 

27. The values of c fi used to form the ratio 

were calculated from the Karman-Schoenherr equation for 
the same Re 0 value for zero heat transfer 


c 


f i 


_ 0.586 _ 

(log 1Q 2 Re Q )(logj 0 2 Re^t-0.868) 


(7) 


The Reynolds number based on momentum thickness was used 
instead of the customary Reynolds number based on the 
distance from the leading edge. A T'—'olds number based 
on a boundary layer parameter was selected because of the 
arbitrariness inherent in calculating an effective leading 
edge for boundary layer measurements on a wind tunnel wall. 

28. The variation of with heat transfer parameter 

(T e - T W )T C is shown in Figure 15 for all of the data of 

the present investigation. Also shown on this figure are 
curves calculated using an extended Donaldson analysis 
reported m reference (u). For those values of Mach number 
for which several data points are available (M « 5.0, 5.8, 
and 6.8), the results indicate that increasing values of 
heat transfer parameter have little effect on the skin 
friction ratio. On the ether hand, Van Driest and Monaghan 
(references (s), (y)) predict an increase in the skin friction 
ratio up to about 10 percent for the highest heat transfer 
case investigated in the present experiments. It should be 
noted, however, that it is implied in these analyses that 
the edge of the laminar sublayer occurs at fixed values of 
u + and y + , and the slope of the turbulent portion of the 
u"**, y* curve remains constant regardless of Mach numlcr or 
heat transfer conditions. As pointed out previously, the 
results given in Table I show that these assumptions are not 
generally valid. This is more evident from Figure 16 in 
which the values of u^ c y^ arc plotted as a function of 

the heat transfer parameter. Also shown are the associated 
theoretical curves derived from reference (u). These figures 
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indicate that Mach number has little influence on the 
variation of c with (T e - T w )/ 7 e * appears, 

however, from these results that caution should be exercised 
in using theories which embody constant values of u^“ * 

for cases of heat transfer. 

28. The variation of with Mach number is shown in 

Figure 17 together with the diiect skin friction measurements 
of Coles (reference (z)) and the deduced skin friction data 
of references (a), (q), (r), (t), (y), (z), and (aa). In 
the preparation of this figure the data for M » 5 and 6.8 
were reduced to the case of zero heat transfer by fitting 
the results of Figure 15 with a straight line using the 
method of least squares. For the other Mach numbers the 
zero heat transfer points were obtained by estimating the 
slope of a straight line 'hrough the data consistent with 
the M *= 5 and 6.8 results. The position of the data 
points on this figure was fixed by plotting the variation of 
c^/Cji against Re^ for each Mach number and selecting 

the values of c^/c fi at a constant value of Re^ of 8,000. 

This value of Ro^ was selected because it represented a 

Reynolds number for which the greatest overlap of data 
existed. It is apparent that the c f /c fi variation with 

Mach number is a smooth continuous curve for a constant 
Reynolds number. Also shown in this figure are the theoreti¬ 
cal predictions of reierences (q) and (u) for a constant 
Reynolds number of 8,000 and zero heat transfer. The close 
fit of the experimental data with the theoretical curves 
indicates that the theories of references (q) and (u) may be 
used to predict with good accuracy the variation of skin 
friction coefficients with Mach number for the zero heat 
transfer case. 


CONCLUDING REMARKS 

30. Detailed investigations of turbulent boundary layer 
velocity and temperature profiles have been made on a nozzle 
wall of tne NOL 12 x 12 cm Hypersonic Tunnel No. 4 at Mach 
numbers of 5.0, 5.8, 6.8, 7.7, and 8.2 for varying rates of 
surface heat transfer. 

31. The experimental results show that turbulent boundary 
layer profiles in hypersonic flow qualitatively resemble 
turbulent boundary layer profiles in incompressible flow in 
many details. 

3L. All velocity profiles measured can be fitted with a 
power law in the outer turbulent portion. 
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33. As in the case of incompressible flow, the turbulent 
portion of the profiles differs in shape from the logarithmic- 
law velocity profile. The discrepancy increases with 
increasing heat transfer. 

34. The dispersing effects of Reynolds number and heat 
transfer parameter exhibited in the u* ( y* velocity profile 
representation for the outer turbulent portion could be 
removed by basing the coordinates on local properties at each 
point in the boundary layer and dividing the distance parame¬ 
ter by the respective shape parameter (II). The data for any 
one Mach number then fall on a single curve. The velocity 
profile curves show a systematic displacement with Mach number. 

35. The velocity profiles in the laminar sublayer are linear. 
The ratio of laminar sublayer thickness to total boundary 
layer thickness (eT L /<y*) u decreases slightly with increas¬ 
ing heat transfer rate but increases considerably with 
increasing Mach number at the same heat transfer rates. 

3G. The u* = y* value at the edge of the laminar sublayer 
for the zero heat transfer case is close to the incompressible 
flow value. However, this n* = y* value is found to increase 
with increasing heat transfer for any given Mach number. 

37. The ratio of the thermal sublayer thickness to the 

total boundary layer thickness ( was found larger 

than (<T L /<n u . The thermal sublayer thickness ratio is 
unaffected by heat transfer but somewhat affected by Reynolds 
number. 

38. The data strongly support the applicability of Reynolds 
analogy for turbulent l^undary layers in hypersonic flow. 

39. Values of skin friction coefficients calculated from 
the heat transfer measurements by using Reynolds analogy 
agree to within 5 percent with those skin friction values 
deduced from the velocity profile slope in the laminar 
sublayer. 

40. For the test range of the heat transfer parameter 

(T 0 - T w )/T c (approximately 0 to 0.5), values of the skin 

friction coefficient are found to bo only slightly affected 
by heat transfer rate and arc in accord with the direct skin 
friction measurements of Coles and other investigators when 
the results are based on a single value of Reynolds number. 

41. The present results extend the range of available skin 
friction and heat transfer data to a Mach number of 8.2. 

42. The detailed measurements of velocity and temperature 
distributions across turbulent boundary layers in hypersonic 
flow presented here enlarge the fund of available data for 
turbulent boundary layers in compressible flows. 
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26.50 8.18 1.000 1.000 
25.30 8.05 1.014 .993 
24.00 7.88 1.045 .985 

21.50 7.30 1.180 .966 
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FIG.4 MACH NUMBER VARIATION ACROSS BOUNDARY LAYER FOR FREE- 
STREAM MACH NUMBERS OF 5 AND 6.8 AND VARIOUS RAI lS OF 
HEAT TRANSFER 
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FIG. 7 VARIATION OF 
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(DATA FROM REFERENCES m, n AND o) 


30 






NAVORD REPORT 3880 








FIG.9 SEMI-LOGARITHMIC REPRESENTATION OF BOUNDARY LAYER 
VELOCITY PROFILES FOR A FREE-STREAM MACH NUMBER 
OF 5.8 AND THREE RATES OF HEAT TRANSFER 
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FIG 10 SEMI-LOGARITHMIC REPRESENTATION OF BOUNDARY LAYER 
VELOCITY PROFILES FOR A FREE-STREAM MACH NUMBER 
OF 6 B AND FOUR RATES OF HEAT TRANSFER 
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SEMI-LOGARITHMIC REPRESENTATION OF BOUNDARY LAYER VELOCITY 
PROFILES FOR FREE-STREAM MACH NUMBERS OF 7.7 AND 8.2 
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FIG. 12 SEMI-LOGARITHMIC REPRESENTATION OF TURBULENT PORTION 
OF BOUNDARY LAYER PROFILES IN u*y*/H COORDINATES FOR 
VARIOUS MACH NUMBERS 
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FIG. 16 VARIATION OF u + L = y + L VALUES 
WITH HEAT TRANSFER PARAMETER 
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FIG. 17 VARIATION OF SKIN FRICTION COEFFICIENT WITH MACH 
NUMBER FOR ZERO HEAT TRANSFER AND A Reo VALUE 
OF 8,000 
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